Vortex-chain phases in layered superconductors 
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Layered superconductors in tilted magnetic field have a very rich spectrum of vortex lattice 
configurations. In the presence of in-plane magnetic field, a small c-axis field penetrates in the 
form of isolated vortex chains. The structure of a single chain is mainly determined by the ratio 
of the London [A] and Josephson [Xj] lengths, a — X/Xj. At large a the chain is composed of 
tilted vortices [tilted chains] and at small a it consists of a crossing array of Josephson vortices and 
pancake stacks [crossing chains]. We studied the chain structures at intermediate a's and found 
two types of behavior. (I) In the range 0.4 < a < 0.5 a c-axis field first penetrates in the form of 
pancake-stack chains located on Josephson vortices. Due to attractive coupling between deformed 
stacks, their density jumps from zero to a finite value. With further increase of the c-axis field the 
chain structure smoothly evolves into modulated tilted vortices and then transforms via a second- 
order phase transition, into the tilted straight vortices. (II) In the range 0.5 < a < 0.65 a c-axis 
field first penetrates in the form of kinks creating kinked tilted vortices. With increasing the c-axis 
field this structure is replaced via a first-order phase transition by the strongly deformed crossing 
chain. This transition is accompanied by a large jump of pancake density. Further evolution of the 
chain structure is similar to the higher anisotropy scenario: it smoothly transforms back into the 
tilted straight vortices. 

PACS numbers: 74.25.Qt, 74.25.Op, 74.20.De 



I. INTRODUCTION 



Layered superconductors have an amazingly rich phase 
diagram in tilted magnetic field. In the presence of the 
in-plane field, pancake vortices generated by the c-axis 
fieldi^ can form a very large number of different lattice 
configurations. Possible structures include the kinked 
lattice^*^, tilted vortex chains^, and crossing lattices 
composed of sublattices of Josephson vortices ( JVs) and 
pancake-vortex stacks4*& In addition to homogeneous 
lattices, phase-separated states may also exist such as 
dense pancake-stack chains sitting on JVs and dilute lat- 
tice in between&2i±fliii or coexisting lattices with differ- 
ent orientations^ Even though considerable progress in 
this field has been made in the last decade, the satisfac- 
tory understanding of the phase diagram has not been 
achieved yet. All these phases probably do realize in dif- 
ferent materials and experimental conditions. However 
finding ground states in tilted field occurs to be a chal- 
lenging theoretical task and it is even more difficult to 
prove experimentally that a particular lattice configura- 
tion does realize somewhere in the phase diagram. 

The main source of richness of lattice structures in 
tilted field is the existence of two very different kinds of 
interactions between pancake vortices in different layers: 
magnetic and Josephson interactions. The key parame- 
ter, which determines the relative strength of these two 
interaction and plays a major role in selecting the lattice 
structures, is the ratio of the two fundamental lengths, 
the in-plane London penetration depth, A = A a fc, and 
Josephson length Xj — 7s, with 7 being the anisotropy 
parameter and s being the interlayer spacing, a = X/Xj. 
One can distinguish two limiting cases which we refer 
to as "extremely anisotropic" case, a < 0.4, and "mod- 



C N;= 

Li 

— — n- 



FIG. 1: Crossing (left) and tilted (right) vortex chains. Up- 
per pictures show three-dimensional views and lower pictures 
show the structures of isolated chains. 



erately anisotropic" case a > 0.7. Note that in our 
terminology even "moderately anisotropic" superconduc- 
tors may have very large anisotropy factor, 7 3> 1. 
Among known atomically layered superconductors, only 
Bi 2 Sr 2 CaCu20 :I ; (BSCCO) and related compounds may 
belong to the "extremely anisotropic" family. Even in 
this compound the parameter a is not smaller than 
~ 0.25 and increases with temperature so that BSCCO 
typically becomes "moderately anisotropic" in the vicin- 
ity of transition temperature. 

In a wide range of the in-plane fields (10-200 G) and 
at very small c-axis fields (up to 1-2 G) the pancake 
stacks in layered superconductors within a wide range 
of anisotropies are arranged in chains, see Fig.^ An iso- 
lated chain is a two-dimensional array of pancake vortices 
oriented perpendicular to the layers. At somewhat higher 
c-axis fields the chains are surrounded by the stripes of 
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regular vortex latticeAiS The internal structure of an 
isolated chain depends on the ratio a and it is relatively 
simple in two limiting cases. At large a the chain is com- 
posed of tilted pancake stacks (tilted chain, right column 
in Fig. ^) and at small a it consists of crossing array of 
Josephson vortices and pancake stacks (crossing chains, 
left column in Fig. ^| . A very nontrivial and intriguing 
problem is how one structure transforms into another 
in the region of intermediate a. We address this prob- 
lem in this paper. We analytically and numerically com- 
puted ground-state configurations in the isolated vortex 
chain and found a surprisingly rich behavior. We found 
two types of phase transitions. The first phase transition 
typically takes place for the intermediate separations be- 
tween pancake stacks a, a = [1 — 2]Xj, and rather wide 
range of the ratio a, 0.4 < a < 0.65. For these a's 
the ground state is given by the crossing chain in a wide 
range of the pancake separations a. However, due to 
attractive coupling between deformed pancake stacks 1 ^, 
the equilibrium separation can not exceed some maxi- 
mum value, which depends on the in-plane field and a 
and it is typically of the order of several Aj. With de- 
creasing the pancake separation o, the crossing chain be- 
comes strongly deformed and smoothly transforms into 
the modulated tilted vortices which then transform via 
a second-order phase transition into the tilted straight 
vortices. We calculated analytically the energies of two 
limiting chain configurations and checked that numerics 
reproduces them. Comparing these energies, we locate 
the transitional region in the phase space, where strongly- 
deformed chains are realized. 

Another phase transition is realized at very small den- 
sities of pancake vortices and only when a exceeds a cer- 
tain critical value rs 0.5 (exact criterion depends on the 
in-plane magnetic field). In this case a small c-axis field 
penetrates in the form of kinks. The kinked vortex lines 
forming tilted chains are composed of pieces of Joseph- 
son vortices separated by kinks. 3,4,5 If the kink energy 
is only slightly smaller then the energy per pancake in a 
straight pancake stack then at very small concentration 
of kinks, typically at a w [20 — 30]Aj, the kinked chains 
are replaced with strongly deformed crossing chains via 
a first-order phase transition. Due to the opposite signs 
of interactions (kinks repel and deformed pancake stacks 
attract each other) this transition is accompanied by a 
very large jump in the pancake density. With further de- 
crease of pancake separation the chain smoothly trans- 
forms back to the tilted chain as it was described in the 
previous paragraph. 

Based on numerical exploration of the chain configura- 
tions, we construct the chain phase diagrams for different 
ratios a. As follows from the above description, there are 
two types of phase diagrams in the region of intermediate 
a's. 

• In the range 0.4 < a < 0.5 a small c-axis field first 
penetrates in the form of pancake-stack chains lo- 
cated on Josephson vortices. Due to attractive cou- 
pling between deformed stacks, their density jumps 



from zero to a finite value. With further increase of 
the c axis field the chain structure first evolves into 
the modulated tilted vortices, which then trans- 
forms via a second-order phase transition, into the 
tilted straight vortices. 

• In the range 0.5 < a < 0.65 a small c axis field 
first penetrates in the form of kinks creating kinked 
tilted vortices. With increasing the c-axis field this 
structure is replaced via a first-order phase tran- 
sition by the chain of pancake stacks, which are 
typically strongly deformed. This transition is ac- 
companied by a large jump of pancake density. Fur- 
ther evolution of the chain structure is identical to 
the smaller a scenario: the structure first trans- 
forms into modulated tilted vortices and then, via 
a second-order phase transition, into tilted straight 
vortices. 

Note that the exact transition between the two types of 
behavior depends also on the in-plane field. As interac- 
tion with the Josephson vortices reduces the energy of 
the pancake stacks, the larger in-plane field favors the 
first scenario. 

Using numerical code developed for studying the chain 
structures, we also investigated stability of an isolated 
crossing configuration of the Josephson vortex and pan- 
cake stack. We found that the crossing configuration 
becomes unstable at a « 0.69. Above this value the 
magnetic coupling is not capable to maintain stable 
configuration and the Josephson vortex tears the stack 
apart. Nevertheless, the obtained stability range oc- 
curs to be significantly broader than it was estimated 
from simple considerations in Ref. 0. The reason is 
that the strongly-deformed crossing configuration sig- 
nificantly modifies the Josephson vortex which reduces 
forces pulling pancakes apart. We found that the crossing 
energy increases smoothly up to instability point. Pertur- 
bative calculation^ gives accurate results for the crossing 
configuration and its energy up to a w 0.35. 

Recently, the vortex chains in BSCCO at small concen- 
trations of pancakes have been studied by the scanning 
Hall probe microscopy by Grigorenko et ali^ They ob- 
served that at very small concentration of the pancakes 
the chains are magnetically homogeneous and separate 
pancake stacks are not resolved. When the external field 
exceeds some critical value of the order of several Oersted, 
crystallites of the pancake stacks are suddenly formed 
along the chain and the flux density in the crystallites ap- 
proximately ten times higher then the flux density in the 
homogeneous chains. Our calculations provide consis- 
tent interpretation for these observations. The magneti- 
cally homogeneous chains are interpreted as kinked/tilted 
chains and formation of crystallites can be attributed to 
the low-density [kinked lines]- [crossing chains] first-order 
phase transition (such interpretation has been proposed 
by Dodgsoni 5 -). The observed large density jump also 
comes out from the theory. 

The evolution of the mixed chain+lattice state with 
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increasing temperature has been studied recently by 
Lorentz microscopy^ It was observed that the pancake 
stacks located in chains smear along the chain direction 
above some field-dependent temperature while the pan- 
cake stacks outside chains still remain well-defined. The 
continuous low-density phase transition from crossing to 
tilted chain found and discussed in this paper provides a 
very natural interpretation for this observation. 

The paper is organized as follows. In Section [H] we re- 
view general expressions for the chain energy. In Section 
lllll we perform analytical calculations of the chain energy 
for the two limiting cases: crossing and tilted chain. In 
the next section llVl comparing energies for the two lim- 
iting configurations, we estimate location of the transi- 
tion region. In scction^we review attractive interaction 
between deformed pancake stacks located on Josephson 
vortices^ and derive general formulas for determination 
of the maximum equilibrium separation between the pan- 
cake stacks and the boundaries of stability region with 
respect to density fluctuations. Section IVT1 contains the 
main results of the paper on numerical exploration of 
the the phase diagram. After discussion of numerical im- 
plementation of the model, we explore stability of the 
isolated crossing configuration. In the next two subsec- 
tions we consider two different phase transitions between 
the tilted and crossing chains and two types of phase di- 
agrams which are realized in the region of intermediate 
parameter a. 



II. ENERGY FUNCTIONAL 

Our calculations are based on the Lawrence-Doniach 
free-energy functional in the London approximation, 
which depends on the in-plane phases </> n (r) and vector- 
potential A(r) 
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are the phase stiffness and the Josephson coupling en- 
ergy with Eq = $0/ (47rA) 2 , A = \ a b and A c are the com- 
ponents of the London penetration depth, 7 = X c /X a b 
is the anisotropy factor, and s is the interlayer period- 
icity. The ratio of the two energy scales determines the 
most important length scale of the problem, the Joseph- 
son length, Aj = 7s = y/ J/Ej. We use the London 
gauge, divA = 0. We mainly address the situation when 
magnetic B inside the superconductor is fixed. The c 
component of the field determines the concentration of 
the pancake vortices n v = B z /Qq inside one layer. The 
in-plane phases <j) n have singularities at the positions of 
pancake vortices R„.i inside the layers, 

[VxVi] z = 2i^<f(r-R ry ). 



Logarithmic divergencies in the vicinity of pancake- 
vortex cores have to be cut at the coherence length £ a ;,. A 
useful approach for superconductors with weak Joseph- 
son coupling is to split the phase and vector-potential 
into the vortex and regular contributions, <fi n — (j)vn + (t>rn 
and A = A.„ + A r . The vortex contributions mini- 
mize the energy for fixed positions of pancake vortices 
at Ej = and give magnetic interaction energy for the 
pancake vortices. One can express this part of energy 
via the vortex coordinates R Mj i. In general, the regular 
contributions may include phases and vector-potentials 
of the Josephson vortices. The total energy naturally 
splits into the regular part F r , the energy of magnetic 
interactions between pancakes Fm, and the Josephson 
energy Fj , which couples the regular and vortex degrees 
of freedom, 



F = F r + F M + Fj 



(3) 
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where the discrete gradient V„</> n is defined as V Tl 
pancakesi 



i n _i_i — (f> n , and Um (R-, n ) is the magnetic interaction between 
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u (r, z) = exp(— z)E\ (r — z) + exp(z)i?i (r + z) . 
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where E\(u) = J°° (exp(—v)/v)dv is the integral expo- 
nent (E\ (u) ?s — 7_e — liiu + u at u <^ 1 with 7_e w 0.5772 
being the Euler constant), r = \/ R 2 + (ns) 2 , and L is a 
cutoff length. The regular and vortex degrees of freedoms 
are coupled only via the Josephson energy in which <f> n 
is the total phase composed of vortex and regular contri- 
butions. 

The discrete layer structure has strongest influence 
on the cores of tilted and Josephson vortices. Inter- 
action contributions to the total energy usually can be 
computed within continuous approximation which de- 
scribes the layered superconductor as a three-dimensional 
anisotropic material. This approximation amounts to re- 
placement of summation in the layer index n in Eqs. 
iJSJ, and 10 with integration in the continuous variable 
z = ns and expansion of cosine in Eq. © • In the contin- 
uous approximation one can derive a very useful general 
result for the energy JIJ (see Ref. Ht)) 

p= $g f ri 3 k (l + Agfc 2 )|S z | 2 + (l + A 2 fc 2 )|S||| 2 

87T J (27r) 3 ( 1 + A 2 fc 2) (1 + A 2 fc2 + A2fc 2 ) 

(8) 

in terms of vorticity S(r) of parametrically defined vortex 
lines Ri(X) 

S(r) = ]T / dX^5(T-B4{X)) 

j 

whose Fourier transform is 

S(k) = £ J dX^eM^ 3 (X)). 

j 
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As we will use this formula only for evaluation of inter- 
action energies between vortex lines, we have to subtract 
from it the logarithmically diverging single- vortex terms. 



In this paper we focus on the structure of an iso- 
lated vortex chain with period a in x direction and pe- 
riod c = Ns in z direction corresponding to the tilt- 
ing angle 9 of vortices with respect to the c axis with 
v = tan# = a/c (see Fig. [TJ. The vertical period is 

fixed by the in-plane field B x , c ~ \J 2$ /(v / 37-Bi). 
For BSCCO 7 ~ 500 and this period is approximately 
equal to 20 layers at B x w 50G. We consider the case 
c <C A and in-plane distances much smaller than A e . 
For this particular problem the general energy given by 
Eqs. (0, 0, |J5J|, and JBJ) can be significantly simpli- 
fied using several approximations: (i) we can neglect 
screening of regular phase and z-axis vector-potential; 

(ii) we consider only one-dimensional displacements of 
pancake rows along the chain, H„ t i — (ai + w„,0, ns); 

(iii) we subtract the energy of straight pancake stacks, 
(B x /$ )eps with £ps « So(lnn + 0.5), allowing us to 
eliminate logarithmically diverging pancake-core contri- 
butions; (iv)we drop the trivial magnetic energy term 
B 2 /8tt which plays no role in selection between different 
chain phases. We will use the chain energy per unit area, 
E = c y (F/V - S 2 /(8tt) - (B z /<f> )e PS ) with V being the 
total system volume. With the above assumptions this 
energy can be represented as 
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where (f> Vfn (r) is the vortex phase variation induced by displacement of pancake rows, «„, from the ideally aligned 
positions 

tan (ir(x + u n )/a) tan (ttx / a) 

Vv,n\x, y;u n ) = arctan — - — — — arctan - — — — , 

tann(7ry/a) tanti (iry / a) 

I 



UMr(u n — u m ,n — m) is the interaction energy between to straight stacks 
the pancake rows per unit length, computed with respect 1 

UMr(u, n) = - > [UM(ma + u, n) — UM(ma, n)] , (10) 
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c y = $o /cB x 3> c is the in-plane distance between chains, 
and L z is the total system length in z direction. 

The energy (jHJ) contains a long-range suppression of 
the Josephson energy accumulated from distances c <C 
y -C c y , that is identical in all chain phases and it is 
convenient to separate this term too. The averaged 2-axis 
phase gradient induced by the Josephson vortex lattice 
is given by 

— — — 2ns ir 2ity 



N Nc v 



Evaluating integral 



Cy/2 ( ( - 2ns 

dy 1 - cos V„0„ - -^—B x y 
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we obtain the long-range Josephson energy, Ej.lr, 



p - P y 
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We will define the local energy, Ei oc , sensitive to the chain 
structure as 



Eioc = E — E. 



j-lr • 



(12) 



This part of energy weakly depends on c y and does not 
diverge for c y — » oo. The result is valid for chains 
separated by distance c y smaller than A c . Similar cal- 
culation can be made for isolated chain separated from 
other chains by distance larger than A c . In this situation 
the integral over y converges on distance A c instead of c y 
leading to result -Ej.lr = EjttX c /sN 2 . We will use this 
result in the analytical calculations of the isolated chain 
energy. 

In calculation of magnetic coupling energy one has to 
take into account periodic conditions for pancake dis- 
placements, Un-\- n — u n . In addition, if one selects z 
axis origin at the center of the Josephson vortex then 
symmetry also requires U— n — —u n . Ground state of the 
vortex chain is determined by the minimum of energy @ 
with respect to pancake displacements and regular phase 
distribution. Two simple limiting cases in Fig. ^ corre- 
spond to (i) u n <C a for the crossing-chain configuration 
and (ii) u n = — <z(l — (n— l/2)/JVj)/2 for the tilted-chain 
configuration. 



III. ANALYTICAL CALCULATIONS OF CHAIN 
ENERGIES 

In this section we compute energy of an isolated vortex 
chain with period a in x direction and period c = N s <C A 
in z direction. In general, there are two approaches to 
compute energies of vortex configurations. Using distri- 
bution of phase and vector potential, the total energy can 
be obtained directly from the Lawrence-Doniach func- 
tional by integration of the local energy. This approach 
is always used in numerical computations. Analytically, 



it is more convenient to calculate the total energy by 
summing up energy of isolated vortices and vortex inter- 
actions. Analytic estimates for energy contributions are 
possible in two limiting cases of weakly deformed cross- 
ing chain and chain consisting of tilted vortices (see Fig. 

Comparison of these energies gives an approximate 
range of parameters where one of these competing con- 
figurations is energetically preferable. In contrast to the 
numerical part, we consider isolated chain separated from 
other chains by distance c y 3> A c . For comparison with 
numerical calculations, it will be necessary to extract the 
local part of energy which is not sensitive to the long- 
range behavior. 



A. Crossing chains 

Energy of crossing-lattices chain per unit area is given 
by sum of pancake-stack (Eps), Josephson- vortex (Ejy) 
and crossing energies (E x ) terms, 



Ecl = Epa + Ejv + E 



(13) 



Both pancake and JV terms are composed of single- 
vortex and interaction contribution, Eps — Ep s + E PS , 
Ejv = E s j V + Ejy. 

We start with evaluating pancake-stack energies. Con- 
tribution from energies of individual stacks to the energy 
per unit length is given by 



E s PS = ^-(ln^ + C v 



(14) 



where C v w 0.497 within the Ginzburg-Landau the- 
ory. Using vorticity of the pancake-vortex chain, 
^( r ) = J2j 3(y)3( x ~~ J a ) corresponding to S z (k) = 
5(k z ) exp (ik x ja), we derive from general formula (JSJ 
the pancake-stacks interaction energy 
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Using relation V . exp (ik x jaj) = ± J2 m 5 ( k x ~ ^ir) and 
integrating over k y we obtain 



Epa = 
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For comparison with the energy of tilted chain, it will be 
convenient to represent this energy in the form 



E PS = — In +1E-U[-— ) 



(16) 



with 



oo , 



C(3)a; 2 /2, x < 0.5 
i _ln 2 + 7B , x >V 
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and C(3) = El 



1.202 



Single- vortex JV energy is given b; 



6,18 



E 



JV 



7c V s 



1.54 



The JV interaction energy can be evaluated in the same 
way as the pancake-stack interaction energy ((TBjl and in 
the limit c <C X it is given by 



E 



JV 



eo /'ttX ^ AirX 
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Therefore, the total energy of the JV array, Ejv = E s j 



Ej V , can be written as 



Ejv = ^(^ + ln C - + Cj V 

7c \ c s 



(17) 



with C JV = 1.54 - ln47r + 7 £ w -0.41. The fist term in 
this formula is the long-range Josephson energy, Ej-lr. 
This term is identical in all chain phases and does not 



determine selection between them. This formula is valid 
for an isolated chain and it is different from the long- 
range Josephson energy of the dense JV lattice, c y < X c , 
given by Eq. i|ll|) . The difference amounts to a simple 
replacement c y — ► ttX c . For comparison with simulations, 
we will use only the local part of energy which is obtained 
from the total energy by subtracting £j_lr. 

Using estimate for the crossing energy of Josephson 
vortex and pancake stack for a = X/"/s < 0. 
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we obtain the contribution from the crossings into the 
energy per unit area 
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Finally, combining results (|rl|) . (fTE|l . ljT7|) and (jT§)l we 

obtain the total energy of crossing chain 



J 



E C l = Ep S 



£o 
a 



ttA 



In- 



4ttA 



\2irX/ 7c 7 



to' 



r 



ln(3.5/a)iV 



(20) 



Subtracting the pancake-stack and long-range Josephson 
energies, we obtain local energy, E l g c L — Ecl ~ Eps — 
E Q irX/jc 2 , 



tilt stiffness of pancake stack and contains magnetic and 
Josephson contribution, sm and ej. Magnetic part has 
been calculated by Clem£ 



E , £ = ? ^ (lnAI + Gn , ) __|_ 5 _ 



(21) 



which we will use for comparison with numerical simula- 
tions. 
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The Josephson contribution to the tilt energy appears 
due to suppression of the Josephson interlayer coupling 
by mismatched pancakes. It can be evaluated as 



B. Tilted chains 

Energy of the tilted chain per unit area also can be 
decomposed into the single-vortex and interaction con- 
tributions 
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TV 
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TV 



E 



TV 



The first term can be estimated analytically only in two 
limiting cases tan# <^ 7 and tan# > 7, where the second 
case corresponds to kinked lines. For the interaction term 
we derive a general formula valid in both limits. 



1. Energy of a single tilted line 

a. Region v = tan 8 < 7 Energy difference between 
tilted and straight pancake stacks, Etv — £ps, determines 



£j — £0 



2 7 2 



hi 



X r , 



Cj 



Numerical constant Cj can be computed exactly by 
matching logarithmic contributions coming from small 
and large distances^ which gives Cj = 3 In 2 — 7^. 
Therefore, the single-vortex energy of the tilted chain 
is given by 
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(22) 

Note that at not very large tilt angles, v < 1, the Joseph- 
son tilt energy is roughly 7 2 smaller than the magnetic 
tilt energy. Two terms become comparable at very large 
tilt angles, v ~ 7. In several papers£*ii the energy of a 
tilted line in anisotropic three-dimensional superconduc- 
tor has been calculated. This calculation is based on the 
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London energy expression JSJ employing elliptic cut off in 
k- integration k 2 + k 2 /j 2 < to treat the vortex core. 
Strictly speaking, this approximation does not describe 
layered superconductors in which £ c < s. Nevertheless, 
in the region v < 7 it leads to the result similar to Eq. 
(|22J) with somewhat different magnetic part and different 
expression under logarithm in the Josephson part. 

b. Region v — tan# > 7 (kinked lines). In the re- 
gion tan# > 7 the vortex lines have the kink structure, 
i.e., they are composed of kinks separated by JV pieces. 
Energy of such line per unit length is given by 

£ki = £ jv + u k /L + e ki (23) 

where e,jv = (eo/7) (l n (V s ) + 1-54) is the energy of 
Josephson vortex, 

u k = se (ln(7*/0 + C k ) (24) 

is the kink energy^, L = st&nO is the separation be- 
tween kinks, and £ki is the kink interaction energy. Nu- 
merical constant C k in the kink energy has been esti- 
mated as —0.17 in Ref. (i. However more accurate nu- 
merical calculations of this paper give somewhat smaller 
value C k w -0.31. 

The interaction energy between neighboring kinks in 
the kinked line decays as 1/L 2 up to L < A c leading to 
relatively large interaction contribution to the total en- 
ergy. The kink interaction energy is computed in details 
in the Appendix lAl In the region L <C A c the dominant 
interaction term in the energy of single line is given by 

I 



Combining all contributions, we derive the follow- 
ing result for the single-vortex energy of a tilted chain, 
E TV — £ki/c, in the limit v > 7 

TPS rps I rps 

£j T y — Upg -t- Ujy 

with Ckv = Cfc — G v w —0.81. Note that the values 
of the numerical constants C v and C'k depend on the 
core structure at small distances r ~ £ from its center, 
which evolves with temperature decrease. However, as 
this structure is exactly the same for the pancake vortex 
and kink, the difference Ck — C v is not sensitive to behav- 
ior at small distances and remains the same down to low 
temperatures. Criterion In {\/a c ) + Ck v = (correspond- 
ing to a c w 0.44) separates the kink and pancake-stack 
penetration regimes of a small c-axis field for large N 
(somewhat larger value a c w 0.5 has been given in Ref. 
6). This critical value increases with decrease of N due to 
the crossing energies in the crossing chains. Penetration 
of the c-axis field in the presence on the in-plane field is 
frequently referred to as a lock-in transition. 

2. Interaction energy of tilted vortices. 

The interaction energy between tilted lines is not influ- 
enced much by the layered structure and it can be com- 
puted within the London approximation. In Appendix iBl 
we derive a general analytical formula for the interaction 
energy of tilted chain, E % TV , 
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1 + VTT. 



(27) 



The first term in this formula represents again the long-range Josephson energy and it is identical to the first term in 
the energy of the Josephson chain Q17jl. At large N the above formula also gives the interaction energy of a kinked 
line for v > 7, because the kinked structure of tilted lines starts to influence their interaction only at very large tilting 
angle v > Nj/2irm& In the region i/«7 Eq. I|27|l simplifies as 
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(28) 



The last term represents change of the Josephson en- 
ergy due to misalignment of pancakes in different stacks. 
Combining this term with the last term in Eq. H22JI. we 
obtain the total Josephson energy loss of tilted chain due 
to pancake misalignment 
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In the limit of kinked lines, v 3> 7, interaction energy 
reduces to 
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Note again that this result is valid until v < N r y/2w, 
where the kink structure of the tilted lines does not influ- 
ence much interaction between them. The limit of larger 
v corresponds to the regime of "kink walls" described in 
Ref . H 



3. Total energy of tilted chains 

Combining the interaction energy (|28[) with the energy 
of individual stacks (|22|l . we obtain the total tilted-chain 
energy in the limits c — Ns <C A, v <C 7 



Etv ~ E' 
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(30) 



This gives the following result for the local part of energy, 
E T v - E PS - EoirX/jc 2 , 



rploc 



C/ (2^) + 2^ (lnAr -°- 95) 



(31) 



In this formula the first term represents the loss of the 
magnetic coupling energy in the tilted chain and the sec- 
ond term represents the Josephson energy loss. 

In the region of kinked lines v ^> 7 (but v <§; Nj/2ir), 
the total chain energy is obtained by combining Eqs. l|2t)|) 
and ijinj) giving 
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(32) 



with Cm 



ln47r - 1/2 + 7B w -2.454. 



IV. LOCATION OF TRANSITIONAL REGIONS 
IN THE PHASE SPACE 

To find out whether the crossing or tilted chain is re- 
alized for given values of the parameters a, c, and a, we 
have to compare the energies of these states. Naively, one 
may think that intersection of the energy curves for the 
two states would correspond to a first-order phase tran- 
sition between these states. However, as we will see from 
numerical simulations, in the region v = tan# < 7 an- 
other scenario is realized. Typically, strongly deformed 
intermediate chain configurations develop in the transi- 
tional region providing a smooth transition between the 
two limiting configurations. Therefore, a simple energy 
comparison gives only an approximate location of the 
transitional region separating the two configurations. 

In the region 1//7 = a/NXj <C 1 comparison of l|20|) 
and (|30|l gives the following criterion for the transitional 
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ln(3.5/a)iV 







(33) 



One can observe that the main competition takes place 
between the loss of the magnetic coupling energy for 
the tilted chain (first term) and strong suppression of 
the Josephson energy by JVs for the crossing chain 
(second term). Solution of this equation provides the 
boundary which can be written in the reduced form 
a = XjN f(N,a) with f(N,a) < 1, i.e., the boundary 
shape in the plane a/Xj-N depends only on the param- 
eter a. 

In region 1 < 1//7 < N/2tt and a > 2irX comparison of 
the energies (|20() and (|32|l gives 
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= (34) 



This equation has a solution only in the kink penetration 
regime, ln(l/a) + Ckv < 0, near the transition between 
the two penetration regimes |ln(l/a) + Ck v \ "C 1 where 
the kink energy is only slightly smaller than the energy 
per pancake of a straight pancake-vortex stack. In con- 
trast to Eq. (|33|l . which gives only an approximate loca- 
tion of the broad transitional region, this equation indeed 
describes a very strong first-order phase transition. 



V. ATTRACTION BETWEEN DEFORMED 
PANCAKE STACKS AND TILTED VORTICES. 
MAXIMUM EQUILIBRIUM SEPARATION 

A peculiar property of the crossing chain is an attrac- 
tive interaction between the deformed pancake stacks at 
large distances^ As a consequence, when the magnetic 
field is tilted from the layer direction, the density of the 
pancake stacks located on the Josephson vortices jumps 
from zero to a finite value. This means the existence of 
a maximum equilibrium separation a m between pancake 
stacks, i.e., chains with a > a m are not realized in equi- 
librium. Note that the tilted vortices also attract each 
other within some range of angles and distances^ meaning 
that tilted chains also have this property in some range 
of parameters. 

A simple analytical formula for the attraction energy 
between the deformed pancake stacks can be derived 
for very anisotropic superconductors Xj 3> A in the 
range A -C R -C XjM- In this limit short-range pan- 
cake displacements u n from the aligned positions in the 
two neighboring stacks produce a dipole-like contribution 
to the interaction energy per unit length between these 
stacks 
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This term has to be combined with the usual re- 
pulsive interaction between straight stacks Uio(R) = 
2e a K (R/X) « 2e a ^irX/2Rexp (-R/X). Minimum of 
the total interaction energy, Uio(R) + SUi(R), gives an 
estimate for the maximum equilibrium separation a„^, 
a rn = Aln(CA 2 / (u 2 )) and this result is valid until a m < 
Xj. Because in BSCCO Xj is only 2 — 3 times larger than 
A, this simple formula is not practical for this compound. 
We will see that a m in BSCCO is usually larger than Xj. 

One can obtain a useful general recipe for determi- 
nation of the maximum equilibrium separation between 
the pancake stacks directly from the chain energy per 
unit area E, without splitting it into the single-vortex 
and interaction parts. We consider situation when the 
in-plane component of the magnetic field is much larger 
than the lower critical field in this direction so that the in- 
plane magnetic induction practically coincides with the 
in-plane external magnetic field. The c component of the 
external field, H z , determines the effective chemical po- 
tential n for density of pancake stacks, u = ^qH z /(Att). 
With fixed chemical potential, the chain thermodynamic 
potential per unit area depends on pancake linear density 



G(n) ~ E(n) — /in. 



(35) 



The equilibrium density is given by the minimum of 
this energy 



dE(n) 
dn 



Substituting this relation back to Eq. (|35|l we obtain 

G(n„) =E{n lx )-E'(n^)n lx . 

If we represent the energy as a function of the stack sep- 
aration a rather than the density, the last result can be 
rewritten in a more compact way 



only up to a certain critical value of separation. Above 
this value the system becomes unstable with respect to 
density fluctuations leading to formation of high-density 
chain clusters. The stability criterion of the chain is given 
by 



d 2 E{n) 
dn 2 



> 0. 



(37) 



As 



d 2 E(n) , d 2 , , , . . 
the stability criterion can also be written as 



d 2 U(a) 
da 2 



> 0. 



(38) 



Therefore, in the dependence of the pancake part of en- 
ergy per unit cell, U, on pancake separation a, the mini- 
mum gives the maximum equilibrium separation and the 
inflection point corresponds to the boundary of instabil- 
ity with respect to cluster formation. 

As an example, we apply the obtained general formula 
(|36[) to the tilted chain at v < 7. Using Eqs. 13Ufl and 
(|17|l . we obtain from Eq. (|36[) a cubic equation for a m 



2| (In N - 0.95) — -jjjr (In N — 0.41) 



m - na = 0, 
. (39) 

with a m — a m /Xj, which can be easily solved numerically 
in general case. At small a and large N an approximate 
solution of this equation gives 

(7TCX \ 
l + _ (In TV _0.41) J , 

i.e., at large N a m approaches a remarkably simple uni- 
versal value 7rA. 



G(a) 



da 



(aE(a)) 



The separation a in this formula corresponds an equilib- 
rium separation only if G(a) is smaller than the energy 
of the Josephson vortex lattice G(oo) = E(oo) = Ejy- 
Therefore the maximum equilibrium separation a m is 
given by the condition 



dU(a) 
do 



= 



(36) 



where the quantity U(a) — a {E(a) — E(oo)) represents 
the pancake part of energy per one stack and the condi- 
tion H36|) implies that a m is determined by the minimum 
of this energy. When the main contribution to the total 
interaction energy is coming from the nearest-neighbor 
interaction, a m coincides with the position of the mini- 
mum in the pair interaction potential. 

In principle, noncquilibrium structures with a > a m 
can be prepared by applying external stretching forces 
at the chain edges. However, the chain can be stretched 



VI. NUMERICAL EXPLORATION OF CHAIN 
STRUCTURES 

A. Numerical implementation of the model 

Our purpose is to calculate the equilibrium distribu- 
tion of the regular phase ip r , n (x, y) and pancake-row dis- 
placements on the basis of the energy To facilitate 
numerical calculations, we introduce reduced variables 

V = v/^Ji x = x/a, v n = u n /a 

and represent the row interaction energy (|10|l in the re- 
duced form 

U Mr (u,n) = ^^V M r(-,n), 
X z a 

where VA/ r (u,n) = VM r (v,n) — Vmi-(0, n) and for 
VAf r (", 1) we derive from Eq. Q 
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(40) 



The reduced energy per unit area, = EXj /e , takes the form 

x r~ a dx r^ 2 



ttN ^ 
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c y /2 



dy 



1 2 

2 (Wr.n) +l-COs(V„ {<t>r,n + 4>v,n)-hy) 



2a 2 N, 



7Tj— ^2 VMr(v n -v m ,n-m) (41) 



717=771 



with = 2ttB x js 2 /&o and a = a/Xj. The reduced local 
energy is defined as 

In particular, from Eqs. (|21|l and (|31|) we obtain the fol- 
lowing results for the reduced local energies of two limit- 
ing configurations in the limit a/N — v/j <C 1 

( hiN +c.jv _ 8 ° 2 for crossing chain 

~ J Naln(3.5/a) , ° 

l ° C ~ U^^j + jlr (In AT- 0.95), for tilted chain. 
' " (43) 



We will also use the excess pancake part of energy defined 
as 

5E = E- E JV (44) 
= (Xj/e )(E-E s PS -E JV ) (45) 

with -Ejv being the JV lattice energy per unit area and 
the excess pancake energy per stack (in units of £o), U = 
a5E. 

The phase distribution ip r _ n and pancake displacements 
v n minimizing the energy functional (|41|l obey the follow- 
ing equations 



A0 n „+sin(V n (<f> r ,n+<l>v,n)-hy)-wx(y n {(j> rtn -x+4> v<n -i)-hy) = (46a) 

oo 

Vy<f>rn(v n a, 0) + —j F M r(v n - v m , n - m) = 0, (46b) 
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where Turiv-, ft) = — VuVMr(^ 5 ft) is the magnetic interaction force between the vortex rows, 

_ , , ttsX ( s\n\\ sX ^ 1 / ^{v-mf + s 2 n 2 

s Mr {v, n)= T exp — cot nv + — } exp — * — . (47) 

a z \ X I a z * — ' v — m I A/a 

m— — oo » 



Note that both Vj^ r (v,n) an d 3~Mr(v,n) are regular at 
v — > 0, because divergency in the first term is compen- 
sated by the to = term in the sum. 

In numerical calculations we assume for simplicity pe- 
riodic boundary condition in y direction with period c y . 
Physically, this corresponds to rectangular arrangement 
of both Josephson vortices and pancake stacks. Even 
though such arrangement does not give the true ground 
state, it does not influence much the structure of an indi- 
vidual chain, that is the focus of this paper. Due to sym- 
metry properties, it is sufficient to find the phase distri- 
bution <fi n (x,y) within domain 0<x<a, 0<y< c y /2, 
1 < n < Ni = N/2 with the following boundary condi- 



tions for the total phase 

<j> n (a + 0,y) = 4> n (+0, y), <t>n(x, -0) = 7T - <j> n (x, +0), 

/ c v \ ( c v \ 2n(n— 1/2) 

tf> n fx, y + 0) = -<p n (x, f - 0J + N , 

<fa(x,y) = -<t>\{-x,y), <j) Nl+1 (x,y) =tt- <j) Nl {x,y). 

Pancake displacements v n have symmetry properties 
Vn+N = Vn and Vi- n = —v n . This allows us to reduce 
infinite to summation in Eq. 146bJ) to summation over 
half of the unit cell 1 < m < Ni , 
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FIG. 2: (a) The dependence of the maximum pancake dis- 
placement in the crossing (defined in the inset) on the ratio 
a = A/A j for different periods a in x direction. Plot shows 
that crossing becomes unstable near a — 0.69. Inset shows 
crossing configuration near instability point, (b) The depen- 
dence of crossing energy on the ratio a. For comparison small 
a calculation is shown in both plots. 



Similar decomposition was made for the magnetic inter- 
action energy in Eq. (|41|) . 

We explored chain phase structures by numerically 
solving Eqs. t4"ffl) with respect to the pancake row dis- 
placement v n and regular phase distribution cf> rtTl for dif- 
ferent values of the parameters a, N = 2Ni, and a. 
Numerical computations were performed on two Linux 
workstations with 2GHz AMD Athlon processors and on 
the nodes of the Argonne computing cluster "Jazz" (350 
nodes, each with 2.4 GHz Pentium Xeon processor). In 
the following sections we review the results of these cal- 
culations. 



B. Stability of crossing configuration 

The fundamental property of the crossing-chain state 
is the structure of the crossing configuration of Joseph- 
son vortex and pancake stack. In-plane currents of the 
Josephson vortex displace the pancakes in the opposite 



directions above and below central layers. Equilibrium 
displacements are the result of a balance between the 
pulling forces of the Josephson vortex, which try to tear 
the pancake stack apart, and magnetic coupling, which 
tries to keep the stack aligned. The structure of the cross- 
ing configuration can be calculated analytically in the 
regime of high anisotropy A <C Xj^i^ leading to the re- 
sult 118|) for the crossing energy and to the maximum dis- 
placement u\ ph 2.2A 2 / [A j In (2Aj/A)]. This calculation 
is based on (i) quadratic approximation for the magnetic 
tilt energy and (ii) assumption that the JV in-plane cur- 
rents are not influenced much by pancake displacements. 
The first assumption breaks down when u\ approaches A 
and the second one breaks down when u\ become compa- 
rable with A j. This means that the both approximations 
break down as A approaches Xj. 

The stability of the crossing configuration has been 
addressed recently by DodgsonSi using the full magnetic 
coupling energy but without taking into account modifi- 
cation of the Josephson vortex by pancake displacements. 
The latter can not be easily computed analytically. This 
calculation suggested that the crossing configuration be- 
comes unstable at a = A/Aj « 1/2.86 w 0.35. This 
estimate seems to be in contradiction with the recent dec- 
oration experiments^ where the isolated pancake stacks 
sitting on the Josephson vortices have been observed in 
the strongly overdoped BSCCO with the ratio a signifi- 
cantly larger than this value. 

To resolve this contradiction and find an accurate sta- 
bility criterion we studied numerically evolution of the 
isolated crossing configuration with increasing ratio a. 
For this purpose we used the code, which calculates the 
chain structure in Fig. ^ for large values of the peri- 
ods a and N = 2Ni . Figure shows the dependence 
of the maximum pancake displacement on the ratio a 
for N — 20 and different values of a. One can see 
that the crossing configuration becomes unstable near 
a 0.69, which significantly exceeds a simple estimate 
111 Ref. 21. The main reason for the extended stabil- 
ity range is that the pancake displacements significantly 
modify the structure of the Josephson vortex. This re- 
duces forces which pull the pancake stacks away and com- 
pensates for reduced magnetic coupling restoring forces 
at large u. Finite-size effect is only noticeable in vicin- 
ity of the instability. The small-a calculation correctly 
predicts the maximum displacement up to a « 0.35 and 
underestimates it at larger a. Figure^ shows the a de- 
pendence of crossing energy e x • It demonstrates a rather 
regular behavior almost up to the instability point. Sim- 
ilar to the maximum displacement, the small-a calcula- 
tion gives an accurate estimate for e x up to a ~ 0.35. For 
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higher a an absolute value of e x exceeds the analytical 
estimate. 



C. Typical phase diagram within the range 
0-4 i$ a, < 0.5: Phase transition from crossing to tilted 
chains with decreasing pancake separation 

At the first stage we studied the evolution of chain 
structures with increasing a = X/Xj for fixed periods 
a and N. As A increases with the temperature and Xj 
is approximately temperature- independent, increase of a 
corresponds to increase of the temperature in real sys- 
tems except that our calculations do not take into ac- 
count thermal fluctuations. For small values of a and 
N we found that the chain structure evolves smoothly. 
An example of such evolution is presented in Fig. [3] for 
N = 14 and a = 1.52Aj. The main plot shows the de- 
pendence of the maximum pancake displacement from 
the straight-stack position u max /a (defined in the inset) 
on the parameter a. We will use this ratio to character- 
ize the chain structure throughout the paper. It changes 
from zero for straight stacks to (1 — l/N)/2 for tilted 
chains. At small a a weakly deformed crossing configu- 
ration is always realized (see the structure at a = 0.25). 
The pancake displacements grow with increasing a and 
the chain evolves into strongly corrugated configurations 
such as configuration for a = 0.4 in Fig. [3] With further 
increase of a, this structure smoothly transforms into 
modulated tilted lines (see the structure for a = 0.46). 
Finally, the last structure transforms via a second-order 
phase transition into the straight tilted lines. For param- 
eters used in Fig. lathis occurs at a = 0.48. The plateau 
in the dependence U max (a) above this value of a corre- 
sponds to the maximum displacement (1 — 1 /N) /2 in the 
tilted chain. 

To compare numerical and analytical calculations we 
plot in Fig. 0Jl the numerically computed a-dependence 
of the local energy (|42|l together with analytical estimates 
(|43|l for N = 10 and a = 3.4Aj. One can see that the 
analytical estimates accurately reproduce numerical re- 
sults for the weakly deformed crossing chain at a < 0.35 
and for the tilted chain at a > 0.5. However, in between 
the numerical study predicts intermediate configurations 
with energies smaller than the energies of the both limit- 
ing configurations. As a result, a naively expected first- 
order phase transition is replaced by a continuous tran- 
sition occurring at significantly larger a. This behavior 
is quite general. We observed it within the broad range 
of the periods, 6 < N < 20, 2 < a/Xj < 5 and the 
ratios 0.4 < a < 0.6. In Fig. 0p we compare location 
of the phase transition into the tilted-chain state in the 
a-a plane for N = 10 with location of the transitional 
region defined by Eq. 1(55)1. The computed transition line 
is always displaced from the estimated boundary in the 
direction of larger a. The transitional region just marks 
the location of the intermediate strongly deformed chain 
configurations. The observed continuous phase transi- 



tion indicates that tilted lines become unstable with de- 
crease of a. It is known that an isolated vortex line 
in anisotropic superconductors is unstable within some 
range of tilt anglesi^^S We have to note that the sta- 
bility criterion of a chain is not identical to the stability 
criterion of an isolated vortex line and requires separate 
study. At large values of a a continuous transition is re- 
placed by a first-order phase transition. However, as it 
was discussed in Sec. due to the attractive interac- 
tion between the pancake stacks, large separations may 
not realize in equilibrium because a is expected to jump 
from infinity to the maximum equilibrium separation a m . 

To find location of a m in the phase space, we calcu- 
lated the evolution of chain structures with changing a 
for fixed a and N. These calculations, of course, repro- 
duce the chain structures and location of the transition 
line of the previous calculation. Following recipe of Sec. 
IVl we calculated the a-dependence of the excess pancake 
energy per unit stack, U(a) = a(E — Ejy), and find 
a m from the minimum location of this energy. Figure 
[S] shows an example of these dependencies for a = 0.5 
and different N. Increasing U(a) at a > a m implies at- 
tractive interaction between stacks at large separations. 
We can see that for a = 0.5 the separation a m weakly 
depends on N and lies in between 2Aj and 3Aj. Note 
that the TV-dependence of the limiting value of U(a) at 
a — ► oo reflects contribution from the crossing energy of 
an isolated pancake stack with the JV array. 

In Fig. we present the chain phase diagrams in the a- 
A-plane for two values of a, 0.4 and 0.5. Solid lines show 
the phase transition into the tilted-chain state (the de- 
pendence a t (N)). One can see that at larger N the tran- 
sition takes place at smaller a. With increasing a this line 
moves higher meaning that the tilted-chain state occupies 
larger area in the phase space. At large a weakly de- 
formed chain configurations are realized, similar to ones 
shown in Fig. [3J for a — 0.25. With decreasing a the 
chain configuration crosses over into a strongly corru- 
gated state. To mark location of this crossover we show in 
the phase diagrams by dashed lines the pancake separa- 
tion at which the maximum pancake displacement u max 
reaches a/4 (plot a r (N)). This crossover can be viewed as 
"reconnection" of pancake stack segments. Dotted lines 
show locations of a m . We see that a m (N) line crosses 
the transition line meaning that at small N a m falls into 
the tilted-chain region and at large N it falls into the 
crossing-chains region. For a — 0.5 we also show the 
analytical estimate for a m for the tilted chain calculated 
from Eq. H39fl . One can see that it agrees very well with 
numerical calculations. For a — 0.5 and N < 8 the tran- 
sition is of the first order. However, the pancake separa- 
tion at the transition lies above a m meaning that it does 
not correspond to the ground state. 

The obtained phase diagrams imply that at small tilt- 
ing angle of the field with respect to the c axis (corre- 
sponding to small a) the tilted chains have lower energy 
than the crossing chains. This is similar to the situation 
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FIG. 3: Main plot shows the dependence of the maximum displacement (defined in the inset) divided by pancake separation 
on the parameters a for TV = 14 and a = 1.52A./. The chain structures are illustrated at marked points. In the configuration 
pictures circles show positions of the pancake vortices and horizontal lines mark locations of the Josephson vortices. One can 
see that the system evolves from weakly deformed chain (a = 0.25) via strongly deformed chain (a = 0.4) to modulated tilted 
chain (a = 0.46). The last structure transforms via a second-order phase transition at a — 0.48 into tilted straight vortices. 
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FIG. 5: The dependence of the excess pancake energy per 
stack in units of eo, U, on the separation between pancakes a 
for q = 0.5 and different TV. The position of minimum of this 
energy corresponds to the maximum equilibrium separation 
a m (marked for TV = 8 plot). 



at higher fields, in the dense lattice, where the crossing- 
lattices state also is expected to transform into the simple 
tilted lattice at small tilting angle of the fields 



FIG. 4: (a) The dependence of the local part of energy of 
the chain configuration on the ratio a for TV = 10 and a — 
3.4Aj. For comparison, we also show analytical estimates 
for the crossing and tilted chain configuration, (b) Phase 
diagram in the plane [pancake separation a]- [ratio a] for the 
fixed separation between JVs, TV = 10. The solid line shows 
computed transition into the tilted chain state. The dashed 
line shows location of the transitional region where the energy 
of the crossing configuration is equal to the energy of the tilted 
configuration. 



D. Typical phase diagram in the range 
0-5 < a < 0.65: Reentrant transition to kinked/tilted 
lines at small concentration of pancake vortices 

At higher values of the ratio a a new qualitative fea- 
ture emerges in the phase diagram. When a exceeds 
the characteristic value, a small c-axis field penetrates 
superconductor in the form of kinks formingkinked vor- 
tex lines (lock-in transition, see e.g., Refs. LJJ0) The 
critical value of a is determined by combination of nu- 
merical constants in the pancake-stack and kink energies 
and it is given by a c = exp C'kv , where the constant Ckv 
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FIG. 6: Chain phase diagrams in the plane a/Xj — N for two 
values of the ratio a, 0.4 and 0.5. The solid line indicate 
phase transition into the tilted-chain phase. The dotted line 
shows the maximum equilibrium separation a m . The dashed 
line shows crossover "reconnection" line at which the maxi- 
mum displacement exceeds a/4 and weakly deformed crossing 
configuration crosses over into the strongly corrugated con- 
figuration. For a = 0.5 we also show a m for the tilted-chain 
phase analytically calculated from Eq. I|39|l . 
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FIG. 7: The dependence of the excess pancake energy per 
stack on pancake separation a for N = 16 and a = 0.6. The 
two branches correspond to the two different starting states at 
large a, crossing chain and kinked lines. One can see that the 
kinked lines have lower energy at very large a, a > 14.81. 
Crossing chain smoothly transforms back into tilted chain 
with decrease of a. The transformation is completed at a 
second-order phase transition point at at = 2.5A./. Chain 
configurations at four points marked by arrows are shown be- 
low. The evolution of chain configuration along this energy 
curve is also illustrated by an animation!— i 



is defined after Eq. 1)26(1 . To find the value of a c , we 
calculated in Appendix 1X1 the energy of tilted chains at 
very large pancake separations a, which allowed us to 
extract the energy of an isolated kink. This calculation 
gives Ckv ~ —0.81 corresponding to a c « 0.44. It is 
important to note that the critical value of a increases 
with decrease of N, due to the increasing contribution 
of the crossing energies to the total energy of crossing 
chain. Very interesting behavior is expected when a is 
only slightly larger than a c . The competing chain states 
have very different interactions: deformed stacks attract 
and kinks repel each other. Moreover, at the same value 
of the c-axis magnetic induction, B z , the kink separation 
is much smaller than the stack separation and the abso- 
lute value of the kink interaction energy is much larger 
than the interaction energy between deformed stacks. As 
a consequence, with increasing B z the total energy of the 
kinked lines rapidly exceeds the total energy of the cross- 
ing chain and the system experiences a first-order phase 
transition into the crossing-chain state. Due to the at- 
tractive interaction between the pancake stacks, the pan- 



cake/kink separation at which the energy curves cross 
does not give the equilibrium separation for the crossing 
chain and the stack separation jumps at the transition to 
a value slightly smaller than the maximum equilibrium 
separation a m . This means that the phase transition is 
accompanied by jump of pancake density and magnetic 
induction, B z . 

This behavior was confirmed by numerical calcula- 
tions. Figure shows a plot of the dependence of the 
pancake energy U on the pancake separation for N = 16 
and a — 0.6. This dependence has two branches, corre- 
sponding to the two different starting states at large a, 
crossing chain and kinked lines. This branches cross at 
a = 14.8A,/ and the the kinked vortex lines have smaller 
energy at larger a. The variations of U at large a oc- 
cur due to the interaction energy and one can see that 
the kink interaction energy is much larger than the in- 
teraction energy of the deformed stack in the crossing 
chain. With further decrease of a, the crossing chain 
smoothly transforms into the tilted chain, as it was de- 
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FIG. 8: The density dependence of the pancake part of the 
thermodynamic potential per unit area, SE — fin (in units 
of eo/Xj) at different chemical potentials fi corresponding to 
different values of the magnetic field strength for N = 14 and 
a — 0.6. Kinks in the curves separate regions of tilted/kinked 
lines (low n) and crossing chains (high n). The equilibrium 
density in units of 1/Aj is given by the global minimum of 
this energy. One can see that at fi ~ —0.152 the system expe- 
riences a first-order phase transition with very large density 
jump. 
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scribed in the previous Section. The second-order phase 
transition for these parameters takes place at a ss 2.5Aj 
somewhat smaller than the maximum equilibrium sepa- 
ration a m ~ 3.44Aj. 

The pancake separation (or pancake density) in the 
chain can not be directly fixed in experiment. Instead, 
the magnetic field strength, H z , fixes the chemical poten- 
tial hh, = ^qH z /{Att), and the equilibrium density 
is determined by the global minimum of the thermody- 
namic potential G(n) — E(n) — finn. To find evolution 
of density with increasing chemical potential, we plot in 
Fig. [8] the density dependencies of the reduced thermo- 
dynamic potential, SE — [in for different [i and represen- 
tative parameters N — 14 and a = 0.6. As the energy 
of isolated stacks is subtracted in SE, the dimensionless 
chemical potential is shifted with respect to its bare value 
and it is related to the magnetic field strength as 

^_ <t> Q (H z -H cl ) 
4ire 

where H c \ is the lower critical field for H||c. We find 
that for selected parameters the transition takes place at 
fi = u t = —0.152. At a < fi t the global minimum falls 
into the region of kinked lines and at \x > m it jumps 
into the region of crossing chain. Note that the density 
value, at which the energy curves cross (kinks in the lines 
near n = 0.1 in Fig. |SJ), is always larger than the lower 
density from which the jump to the high-density state 
takes place. At the transition, the density jumps almost 
ten times, from 0.037/Aj to 0.315/Aj. 

The numerically obtained phase diagram in the N-a 
plane for a = 0.6 is shown in the upper panel of Fig. 
E3 The plots in the lower left panel, the thermodynamic 
potential at the transition point and the maximum dis- 
placement versus a/A,/, illustrate definitions of different 



FIG. 9: The upper plot shows the phase diagram in the N-a 
plane for a — 0.6. The left plot in the lower plat illustrates 
meanings of the phase lines using the plot of the thermo- 
dynamic potential SE — un vs a/X.j at the transition point 
(main plot) and the maximum displacement u divided by lat- 
tice spacing a vs a/Xj (inset). The lines 1 and 1' correspond 
to the two limiting pancake separations at the transition point 
between which the jump occurs. The line 2 indicates crossing 
of the energy curves for the kinked and crossing chain. The 
line 2' shows the position of a continuous transition into the 
tilted chain (illustrated by the inset in the lower panel). Dot- 
ted line slightly above l'-line shows the position of the maxi- 
mum equilibrium separation a m . We also show the crossover 
line a/Xj — N above which well-defined kinks appear. The 
right plot in the lower part shows blowup of the phase diagram 
above the 2' line. 

lines in the phase diagram. The lines 1 and 1' show the 
limiting pancake separations at the first-order transition 
between which the jump takes place. When the chem- 
ical potential is fixed by external conditions, the area 
between these lines is bypassed in equilibrium. At large 
N the jump takes place from the kinked-lines state into 
the strongly corrugated configuration. This configuration 
transforms into the tilted chain with further decrease of 
a via continuous transition shown by the line 2'. Below 
N = 14 only tilted chains realize, but the density jump 
still exists. The upper separation grows approximately 
proportional to N while the lower separation slowly de- 
creases with N and lies slightly below the maximum equi- 
librium separation a m shown by dotted line. This means 
that the relative density jump increases with N . The line 
2 shows the position of the crossing of the energy curves 
for the two states. 

It is also instructive to examine the phase diagram for 
fixed c-axis period N, in the a-a plane. As a increases 
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FIG. 10: Phase diagram in a-a plane for N = 14. Solid nose- 
shaped line shows the phase transition into the tilted-chain 
state. Dashed line shows location of the maximum equilib- 
rium separation, which terminates at some point. We also 
show by dotted line location of the stability boundaries ob- 
tained from inflection points of the dependencies U(a) (see 
criterion 
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FIG. 11: Left panel: The dependencies of the excess pancake 
energy per stack U on separation a for N = 14 and different 
values of a. With increasing a the first-order transition from 
tilted to crossing chain vanishes and the crossing chain do not 
realize at all. Also the dependence U[a) becomes monotonic 
which corresponds to termination of a m (a) line in Fig. I1UI 
Right panel: The dependence of the relative maximum dis- 
placement u on separation a for N — 14 and a — 0.63. The 
maximum value of u, corresponding to tilted chains, is given 
by (1 — l/14)/2 w 0.464. One can see that there is a range of 
a where the chains become only slightly modulated. 



with temperature, this diagram to some extent describes 
the temperature evolution of the chain structure. Fig- 
ure E3 shows such diagram for TV = 14. One can see 
that the phase transition line is reentrant, within some 
range of a the tilted chains are realized both for small 
and large a. Above some critical value of a (0.66 for 
N = 14) only straight tilted lines exist in whole range of 
a. This maximum value increases with N. For somewhat 
smaller values of a, there is narrow range of a where the 
chains become slightly modulated. This is illustrated in 
the plot of the maximum displacement u for a = 0.63 
in the right panel of Fig. 1111 The maximum equilibrium 
separation line a m (a) terminates at certain value of a, 
above which the dependence U(a) becomes monotonic. 
We also show in this diagram the stability boundaries 
extracted from the inflection points in the U (a) depen- 
dencies (see criterion (|38ll and related discussion). Ori- 
gin of such phase diagram can be better understood by 
studying the dependencies U (a) at different a shown in 
the left panel of Fig. ^2 One can see that the first-order 
transition vanishes above certain value of a and the de- 
pendence U(a) becomes monotonic close to this value of 
a. On the other hand, even after becoming monotonic, 
this dependence still has two inflection points bounding 
unstable region in some range of a. Existence of such un- 
stable region means that the pancake-density jump with 
increasing the c-axis external field persists in the region 
where only tilted chains exist. This jump is closely re- 
lated (but not identical) with the jump of the tilt angle of 
the vortex line at the lower critical field with increasing 
tilt angle of the external field i 27 i 28 

The vortex chains in BSCCO at small concentrations 
of pancakes have been studied by the scanning Hall probe 
microscopy m Ref. B It was found that at very small 
concentration of pancakes the chains are magnetically ho- 



mogeneous and separate pancake stacks are not resolved. 
When the external field exceeds certain critical value, 
crystallites of the pancake stacks are suddenly formed 
along the chains and the flux density in crystallites ap- 
proximately ten times higher then the flux density in ho- 
mogeneous chains. The [kinked lines]- [crossing chains] 
first-order phase transition provides a very natural ex- 
planation for this observations. 

VII. CONCLUSIONS 

In conclusion, we investigated numerically and analyt- 
ically the phase diagram of an isolated vortex chain in 
layered superconductors. In the region where Josephson 
and magnetic coupling are approximately equal, we found 
a very rich behavior. The crossing chains typically trans- 
form into tilted chains with decreasing pancake separa- 
tion via formation of intermediate strongly deformed con- 
figurations and a continuous phase transition. When the 
relative strength of the Josephson coupling exceeds some 
typical value, the phase diagram becomes reentrant. At a 
very small c-axis field, tilted chains are realized in which 
the vortex lines have the kinked structure. With increas- 
ing c-axis field these low-density tilted chains transform 
via a first-order phase transition into strongly-deformed 
crossing chains. This transition is accompanied by a large 
jump of the pancake-vortex density. With further in- 
crease of the field these crossing chains transforms back 
into the tilted chains via a second-order transition. 

An important feature of real BSCCO which is not 
taken into account in this paper is the thermal vor- 
tex fluctuations. We expect that fluctuations will not 
change qualitatively the described behavior, especially 
the strong first-order phase transition, but may signifi- 
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cantly change locations of the phase transitions in the 
phase space. 

Finally, we briefly overview the relevant field scales. 
In the range of studied JV separations, N, from 10 
to 20, and for 7 = 500, the in-plane field, B x = 
2<I>o/(v / 37s 2 iV 2 ), varies in the range from 200 to 50 Gauss 
and the in-plane separation between Josephson vortices, 
c y = (\^3/2)-fsN, varies from 6.75 to 13.5 fim. The 
typical density jump for N = 16 and a = 6 in Fig. 
corresponds to jump if the average c-axis induction, 
B z = &o/ac y , from 0.12 to 1.2 Gauss. Such jump is usu- 
ally difficult to notice in the global magnetization mea- 
surements. Therefore, a rich spectrum of transformations 
discussed in this paper takes place in the range of very 
small c-axis magnetic induction, not exceeding few gauss. 
Taking a typical value A for BSCCO at 80K, as 0.4 /xm, 
we estimate for the same parameters that the maximum 
c-axis field in the chain center, B z q = $o/aA, jumps from 
3.3 Gauss to 33 Gauss. 
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APPENDIX A: KINK INTERACTION ENERGY 
OF A SINGLE LINE 

We compute the kink interaction energy within London 
approximation. As the main contribution to this energy 
comes from the regions away from the JV and kink cores, 
one can expect that the London approach gives a very 
good approximation of the interaction energy. Shape of 
the kinked line is given by 

R(X) = (X,0,u(X)), 

u(X) = ns, (n - l/2)L < X < (n + 1/2)L. 

From the general formula © we obtain the total energy 
of the kinked line in London approximation 



8nL 7 



d 3 k 

(2^)3 



dX / dX 



1 + A 2 fc 2 + A 2 fc 2 



{l + \ 2 k 2 ) ( 



exp(ik x (X - X') + ik z (u - u')) 



(Al) 



To separate the kink interaction energy one has subtract from this expression the energies of Josephson vortices, £j V , 
and kinks, e£ ink , in London approximation. Integration over X , X' , and k x leads to the following expression for the 
total kink contribution to energy, = e^i — £jy 



with 
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g = LJ\ C ~ + 1 ~ 2 k 2 z +k 2 y , 



gi = LyJX- 2 + k 2 + k 2 . 

To separate kink interaction we have to subtract the contribution coming from isolated kinks, i.e., the term which 
behaves as 1/L at L — > 00. 
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This gives the following result for the kink interaction energy 
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which has to be evaluated in the limit L > 7s. The main contribution is coming from the first term in square brackets. 
If we keep only this term, than the kink interaction energy can be reduced to the following form 



with 



J(r) = 



dpydp z 



p =p v 



,2 
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\ /l +P 2 exp {ry/l+p^j - 1 



r 



In the practically interesting case r 
integral J(r) can be evaluated as 



L/X c < 1 the 
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4tt 
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giving the main result for the kink interaction energy 



The second term in square brackets of Eq. 1A2|) repre- 
sents magnetic coupling contribution to the kink interac- 
tions. We calculated this contribution in the two limiting 

cases 



-(2) 



Sj^, for L < A 
^fexp(-L/A), for L< A 



As we can see, it does give a very small contribution to 
the total kink interaction energy. 



APPENDIX B: DERIVATION OF THE 
INTERACTION ENERGY OF TILTED 
VORTICES. 



The interaction potential between two straight tilted 
vortices per unit length along c-axis separated by dis- 
tance R in the tilt direction (x axis) is given by 



$ 2 f dk x dk y v 2 + (1 + \ 2 c {kl + v 2 k 2 x ))/{\ + \ 2 (k\ + v 2 k 2 x )) 



1 + \ 2 V 2 k x 

I 



\lk\ 



exp(ik x R) 



This formula works also in the regime of kinked vortex Interaction energy of the chain per unit area E l TV can be 
lines v > 7. The kink structure of the lines starts to represented as 
influence interaction between them when kink separa- 
tion L exceeds c z j/2ir corresponding to the condition 
v > N-f/2ir. Integrating over k y we obtain Ui{R) = 
J (dk/2n) cos (kR) U l (k) with 
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The last term is small and will be dropped in further 
calculations. Collecting terms, we finally obtain for the 
total interaction energy J2ZJ- 



FIG. 12: Left panel shows plot of the pancake part of the 
chain energy per unit area in units eo/Aj, 5E, vs pancake 
density n at very small n. Slope of this energy at n — > 
determines the energy of an isolated kink and the numerical 
constant Ckv Calculations were made for two sets of param- 
eters, (N = 10, a = 0.52) and (N = 14, a = 0.6), and for 
two system sizes for each set. Right picture illustrates the 
numerically calculated kink structure. Arrows show in-plane 
currents and gray level plot codes distribution of the cosine 
of phase difference between the neighboring layers (dark area 
below corresponds to the Josephson vortex.) 



In the limit a -C 2ir^v 2 \ 2 + A 2 the first term, j — 1, can 
be evaluated as 



pi ~ £ ° 

£j tvx 
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(^ 2 A 2 + Ag) 



y/^ 2 + 7 2 L I 47T^/ h> 2 \ 2 + A 2 
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APPENDIX C: ENERGY OF ISOLATED KINK 



To find the energy of an isolated kink, we calculated 
energy of tilted lines in the regime when kink separation 
L = a/N significantly exceeds the Josephson length. Nu- 
merically, this is a challenging task because the kink in- 
teraction energy decays slowly with increasing L meaning 
that one has to go to very large values of a. To main- 
tain sufficient accuracy, one has to use large number of 
grid points in x direction. As follows from Eq. I|32(l . the 
pancake part of energy vanishes linearly at small kink 
concentrations n. In reduced units, we define this energy 
as 8E = (Xj/e ){E T v ~ E s pv - E JV ) and from Eq. ^ 
we have 



5E sa n ( In - + C kv 

a 



Using integral 
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cik 



1 Va 2 - b 2 + a 

o Vl + a 2 k 2 (l + b 2 ^ 2 ) a>b Va 2 - b 2 n b 

other two terms are calculated in the limit c <C 27rA as 
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a 1 + Vl + v 2 
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c \ 2 I 7 _ 



Plots of this energy are shown in Fig. ^] for two sets of 
parameters, (N = 10, a = 0.52) and (N = 14, a = 0.6). 
From linear fits at small n we obtain estimates Ckv ~ 
-0.192 - ln(l/0.52) » 0.846 for the first set and C kv w 
-0.296 - ln(l/0.6) w 0.807 for the second set. If we 
use the last constant, corresponding to the larger system 
size, we obtain the estimate f=a —0.31 for the constant 
in the kink energy within Ginzburg-Landau theory 124|) . 
This is somewhat smaller than the value —0.17 reported 
m Ref . 0. The difference is most probably due to finite- 
size effects. 
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